TDDFT WITH SKYRME FORCES: EFFECT OF TIME-ODD 
DENSITIES ON ELECTRIC GIANT RESONANCES 



V.O. Nesterenko\ W. Kleinig J. Kvasil ^, P. Vesely ^, and P.-G. Reinhard ^ 

^ Bogoliuhov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, 
Duhna, Moscow region, 141980, Russia 
^ Technische Univirsitat Dresden, Inst, fiir Analysis, D-01062, Dresden, Germany 
^ Institute of Particle and Nuclear Physics, 
Charles University, CZ-18000 Praha 8, Czech Republic and 
'^Institut fiir Theoretische Physik II, 
Universitdt Erlangen, D-91058, Erlangen, Germany 

Abstract 

Time-odd densities and their effect on electric giant resonances are investigated within the self- 
consistent separable random-phase-approximation (SRPA) model for various Skyrme forces (SkT6, 
SkO, SkM*, SIII, SGII, SLy4, SLy6, SkI3). Time-odd densities restore Galilean invariance of the 
Skyrme functional, violated by the effective-mass and spin-orbital terms. In even-even nuclei these 
densities do not contribute to the ground state but can affect the dynamics. As a particular case, 
we explore the role of the current density in description of isovector El and isoscalar E2 giant 
resonances in a chain of Nd spherical and deformed isotopes with A=134-158. Relation of the 
current to the effective masses and relevant parameters of the Skyrme functional is analyzed. It is 
shown that current contribution to El and E2 resonances is generally essential and fully determined 
by the values and signs of the isovector and isoscalar effective-mass parameters of the force. The 
contribution is the same for all the isotope chain, i.e. for both standard and exotic nuclei. 
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I. INTRODUCTION 

Skyrme forces are widely used for exploration of the ground state and excitations 

in atomic nuclei. In particular, last years big effort was done in application of these forces 
to exotic nuclei, nuclear matter and astrophysical problems (for recent reviews see js, 4|). 
However, in spite of the impressive progress, Skyrme forces still suffer some serious troubles. 
For example, there are too many Skyrme parametrizations and the forces do require a 
unification. Further, being successful in description of the nuclear ground states, Skyrme 
forces often fail in a comprehensive treatment of nuclear dynamics (e.g. no one force can 
describe equally well both isoscalar and isovector modes, see jsl for recent discussion). So, 
the further upgrade of the forces is necessary and nuclear dynamics should play here an 
important role. 

In this connection, we propose here the analysis of time-odd densities in the Skyrme 
functional and their role in description of electric giant resonances (GR). Time-odd densities 
(current j, spin s, and vector kinetic-energy T) are known to restore Galilean invariance of 
the Sky.„.e functional, vWated by velooty-dependent tin.e-even densities (kinat.c-anergy . 
and spin-orbital ^3) [g, 17[ . Therefore time-odd densities are related to important features of 
the functional (effective masses and spin-orbital interaction) and might play an essential role 
in its upgrade. In this paper, we discuss time-odd densities in general and, as a particular 
case, scrutinize the influence of the current density j on the isovector E1(T=1) and isoscalar 
E2(T=0) GR in a chain of Nd isotopes with A=134-158. As compared with the previous 
studies p, luj, much more nuclei and forces are involved to make the exploration indeed 
systematic. We will establish important relations of the GR properties to specific parameters 
and terms of Skyrme forces and, on these grounds, propose a tentative classification of the 
forces. 

The calculations are performed within the self-consistent random-phase-approximation 
SRPA) method with factorized Skyrme forces j^, sl, y, 1^, 12 1. SRPA covers both spherical 
9] and deformed j^, isl nuclei. Factorization minimizes the computational effort minimized 
and makes possible the systematic exploration. 
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II. TIME-ODD DENSITIES 



A. General properties 

The Skyrme functional [l[ in the form 



(1) 



includes kinetic, Coulomb, pairing and Skyrme terms respectively. Expressions for the first 
three terms are done elsewhere 
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12, 



13|. The Skyrme part reads 
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This part involves time-even (nucleon pg, kinetic-energy r^, and spin-orbital Q^g) and time- 
odd (current jg, spin Sg, and vector kinetic-energy Tg) densities, where q denotes protons 
and neutrons. The total densities (like j = jp + jn) are given in ([2]) without the index. 

Both time-even and time-odd densities follow from the original Skyrme forces it]. Time- 
even densities contribute to both ground state and excitations. Time-odd densities affect 
the excitations as well. They also influence g.s. of odd and odd-odd nuclei but are irrelevant 
for g.s. of spin-saturated even-even nuclei. 

As was mentioned above, time-odd densities restore Galilean invariance of the functional, 
violated by velocity-dependent time-even densities and ^g Just for this reason, 

time-odd densities enter the functional only in the specific combinations with their time- 
even counterparts. Hence, time-odd densities do not lead to any new parameters. This is 
especially useful for even-even nuclei where Skyrme parameters, being fixed for the g.s. with 
implementation of the time-even densities only, can be further applied to nuclear dynamics 
involving time-odd densities as well. 

Explicit expressions for time-even and time-odd densities via single-particle wave func- 



tions can be found elsewhere 



13| . It is more instructive to write down these 
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densities in terms of the basic scalar pg and vector Sq densities 

Pg{f,r) = ^pifa^a'), s^ify) = ^ pg{faya'){a'\a^\a) (3) 

cr aa' 

which form a general one-body density matrix {sl 

Pqiraya') = hpq{ry)6^y +'^{a\a^\(r')s''g{f,r)] . (4) 
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V 

Here a,, is a Pauli matrix and cr = ±1 . Then the densities in (12]) read as 



Pq{r) = pg{f,r), Sg{r) = Sg{f,r), (5) 

^[(V - V')Pg{f,f^)],=,,, 3f (f) = i 

rg{r) = [V • Vpgir, f')],=r^, f,(f) = [V ■ V)sg{r, r')],=r^ , (7) 



li^ = 7^[(V - V')p,(r,r-^)].-=,,, 3f (f) = -[(V, - V;).^(r,r-)],-.=,,, (6) 



i.e. are reduced to the basic densities ([5]), their momenta and kinetic energies ([7j). This 
presentation clarifies physical sense of the densities in terms of hydrodynamics. Besides, 
the presentation ©-(ITj) can be treated as some kind of the gradient expansion up to the 
second derivatives of the basic densities. Such expansion is obviously useful for non-uniform 
densities, which is the case for Sg and pg at the nuclear boundary 3|. 

B. Current density 

Between time-odd densities, the current is most important for electric GR. Indeed, the 
spin density Sg is mainly relevant for magnetic modes and does not influence electric GR 



The density Tg can also be neglected as it supplements in ([2]) the term which is 
commonly omitted in most of the Skyrme forces. 

Contribution of jg to the residual interaction is driven by the variation j^, 0] 

= 2[-&i + b[6g,g>W - r) (8) 

which is fully determined by the terms ~ bi,b[ in the functional ([2]). In these terms, the 
current density adjoins the values pr and PgTg responsible for the effective mass. Hence, one 
may expect the correlation between effective masses and jg in GR dynamics. 

To inspect this correlation, it is convenient to express the terms ~ hi, b[ through isoscalar 
and isovector densities po = Pn + Pp and pi = Pn~ Pp (the same for r and j). Then the sum 
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Skyrme force 

FIG. 1: Isoscalar Bq and isovector Bi force parameters as well as related isoscalar effective masses 
m* and isovector enhancement factors n for the Skyrme forces SkT6 (1), SkO (2), SkM* (3), SGII 
(4), SIII (5), SLy4 (6), SLy6 (7) and SkI3 (8). The dotted line at the bottom plot marks zero 
value. 



of the terms with hi and h\ is transformed to the form 

Bo{poTq - il)- Bi{piTi - Ji) 



(9) 



where isoscalar and isovector contributions are decoupled and Bq = bi —b[/2 and Bi = b[/2 
^. For the symmetric nuclear matter, the isoscalar effective mass m*/m and the sum-rule 
enhancement factor k, = {m\/m)^^ — 1 (where m\/m is the isovector effective mass) can be 
expressed via the new parameters as 



,m*, 1 m ^ 2m ,^ 

(-)- = l + 2-p5o, n=-piB. 



(10) 



where m is the bare nucleon mass and p is the density of the symmetric nuclear matter. The 
relations (fTOl) are illustrated in Fig. 1 for the set of Skyrme forces used in this paper. We see 
that there indeed take place correlations Bq ^ m* /m and Bi ^ n. In the Section 4 we will 
show that jq-impact to E2(T=0) and E1(T=1) GR is also determined by the parameters Bq 
and Bi. This will justify the relation between jg and effective masses. 



The velocity- and spin-dependent densities in the Skyrme functional should influence 
energy- weighted sum rules (EWSR) Q, Q. For electric modes, contributions of spin- 
dependent densities 3, s and T to EWSR can be neglected (though for Skyrme forces with 



a low m*, like SkI3, influence of the spin-orbital density can be noticeable and reach 7-8% 
jioll). Instead the contributions of r and j are essential. The former alone results in the 
effective mass. However, for isoscalar modes the r and j contributions fully compensate 
each other 15|. So, if we take j into account, then EWSR(T=0) again acquires the bare 
mass m and becomes 



EWSR{T = 0, A > 1) = ^^A(2A + 1)''A < r^^-^ 



iirm 



>A 



(11) 



Such compensation is not complete for isovector modes [15| and EWSR(T=1,A > 1) retains 
the effective mass. In particular. 



EWSR{T = 1,A = r 



8nml A 



(12) 



III. CALCULATION SCHEME 



SkO 



The study was performed with t he representative set __of 8 Skyrme forces: SkT6 [1 



18|, SIII [19|, SGII [20(, SkM* [2l(, SLy4 



2^, SLy6 



22 , SkI3 



23|. The isovector and 



isoscalar parameters and characteristics o: 



Fig. 1. We used the SRPA approach 



these forces relevant for our study are exhibited in 



12| in the approximation of 5 input operators 



of a different radial dependence. Four operators were chosen following prescription [9|]. For 
deformed nuclei, the operator r^~^'^Yx+2,^J. was added to take into account the multipole mixing 
of excitations with the same projection fi'^ and parity tt. Direct and exchange Coulomb 
contributions to the residual interaction were taken into account. As was shown in our 
previous studies , 111! , the proper choice of the input operators results in sensitivity of the 
separable residual interaction to both surface and interior dynamics. Hence high accuracy of 
the calculations was achieved already for a few separable terms. Factorization of the residual 
interaction results in drastic simpliflcation of RPA calculations, which allows systematic 
calculations even for heavy spherical and deformed nuclei. It should be emphasized that 
factorization in SRPA is fully self-consistent and does not lead to additional parameters. 
The giant resonances were computed as energy-weighted strength functions 



(13) 



smoothed by the Lorentz function C,{ijj—ujy) = A/[27r((£'— £',^)^ + (A/2)^)] with the averaging 
parameter A. In all the calculations A= 2 Mev which was found to be optimal for the 



comparison with experiment and simulation of various physical smoothing factors. Further, 
Mxn^ is the matrix element of EXfi transition from the ground state to the RPA state >, 
is the RPA eigen-energy. We directly compute the strength function (IT^ and hence 
fully avoid determination of the numerous RPA eigen-states. This additionally reduces the 
computation time. SRPA calculations of GR in one nucleus at a familiar laptop take less 
than one hour, as compared with weeks while using full RPA methods 24 1. 

The E1(T=1) and E2(T=0) resonances were computed with the proton and neutron 
effective charges e^-^-^ = N/A, e^J^ = —Z/A and e^^-^ = e'^J^ = 1, respectively. A large 
configuration space was used. E.g. in ^^°Nd we took into account 247 proton and 307 
neutron single-particle levels ranging from the bottom of the potential well up to ~ +20 MeV. 
This results in 13000 and 22000 2-quasiparticle configurations for E1(T=1) and E2(T=0) 
excitations, respectively, with the energies up to 65 MeV. Then, EWSR ffTTl) and f|T2|) are 
exhausted by 93-99% and 95-97%, depending on the force. 



IV. RESULTS AND DISCUSSION 



Results of the SRPA calculations are presented in Figs. 2-7. In Figs. 2 and 3 the 
resonances E1(T=1) and E2(T=0) in ^^°Nd are exhibited for eight Skyrme forces. It is seen 
that all the forces give in general acceptable agreement with the experiment for E1(T=1) 
resonance, which confirms ability of SRPA to describe isovector modes. For E2(T=0) we 
have no experimental data and so give the empirical estimation. It is seen that all the forces, 
for exception of SkT6 and SkO, overestimate the E2(T=0) energy and the lower the isoscalar 
effective mass m* of the force, the more the overestimation. Figs. 2 and 3 demonstrate that 
no one force can describe simultaneously E1(T=1) and E2(T=0). This is especially the case 
for the forces with m* /m < 0.8. 

Figures 2 and 3 compare results obtained with and without j-contribution to the residual 
interaction. The careful inspection of the figures shows that j-contribution always improves 
description of E2(T=0) but not of E1(T=1). Moreover, the current always results in the 
down shift for E2(T=0) but may cause various shifts for E1(T=1). By comparing Figs. 2 
and 3 with Fig. 1, one can see that j-effects in E2(T=0) and E1(T=1) are fully determined 
by parameters Bq and Bi, respectively. The isoscalar parameter Bq is positive for all the 
forces and steadily increases from SkT6 to SkI3. Hence we have for E2(T=0) the shift of 
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FIG. 2: E1(T=1) giant resonance in i^°Nd, calculated with the Skyrme forces SkT6, SkO, SkM*, 
SGII, SIII, SLy4, SLy6 and SkI3 for the cases with (solid curve) and without (dotted curve) 
contribution of the time-odd current. The strength is smoothed by the Lorentz weight with the 



averaging A=2 MeV. The experimental data 



the same sign for all the forces and magnitude of the shift rises from SkT6 to SkI3. Instead, 
the isovector parameter Bi is zero and small for SkT6 and SkO, positive for SGII, SIII 
and SkM* and negative ifor SkT4, SkT6 and SkI3. Hence we have negligible (SkT6), up 
(SkM*, SGII, SIII) or down (SkT4, SkT6, SkI3) shifts of E1(T=1). Magnitude of the shifts 
obviously correlates with the absolute values of Bi. 

The above analysis allows to separate the Skyrme forces into 3 groups, depending on their 
Bq and Bi parameters. In the first group (SkT6, SkO) these parameters are zero or small, 
the effective masses are equal or close to the bare nucleon mass, and the role of the current 
density is negligible. These forces usually well describe E2(T=0) but fail for E1(T=1). 
The second group (SGII, SIII, SkM*) has larger Bq and Bi and therefore smaller effective 
masses. What is important, Bi here is positive. These forces are usually fitted by GR 
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FIG. 3: The same as in Fig. 1 for E2(T=0) giant resonance. The empirical estimation for the 
resonance energy, E = 62^-1/3 MeV, is marked the arrows. 

as well (without taking into account time-odd densities). Hence they provide a reasonable 
description of GR. In this case, adding j-contribution can even worsen the description, as 
we see in Fig. 2. These forces violate Galilean invariance and so can hardly be considered 
as a robust basis for the further upgrade. The third group of the forces (SLy4, Sly6, SkI3) 
is characterized by low effective masses and, what is important, by the negative value of Bi. 
These forces (at least SLy4 and Sly6) well describe E1(T=1) but essentially overestimate 
the E2(T=0) energy. No one of the three groups is perfect. More systematic work with 
Skyrme forces is necessary. 

Now we turn to systematic exploration of E1(T=1) resonance in the chain of Nd isotopes. 
The calculations were performed with the force SLy6 which, following Fig. 2 and our previous 
studies , promises one of the best descriptions of E1(T=1). Fig. 4 represents E1(T=1) 

in the isotopes where we have the photoabsorption experimental data. The isotopes run from 
the spherical semi-magic ^^^Nd to the axially deformed ^^°Nd discussed above. The /i = 



9 




15 20 25 



10 15 20 25 30 



E [MeV] 



FIG. 4: The dipole giant resonance in Nd isotopes with A=142-150, calculated with the Skyrme 
forces SLy6. The left panels exhibit the full strength (solid curve) calculated with the current 
contribution. The /x = (left small peak) and /U = 1 (right big peak) branches of the resonance 



27l | are depicted by triangles. The 



are exhibited by dotted curves. The experimental data 
right panels compare the full dipole strengths calculated with (solid curve) and without (dotted 
curve) current contribution. 



and /i = 1 brunches of the resonance are exhibited to illustrate the magnitude and kind 
(prolate) of the deformation. The figure demonstrates a remarkable agreement with the 
experiment for all the nuclei, maybe for exception of ^^^Nd. Following right plots of the 
figure, the current contribution is the same for all the isotopes. 

In Fig. 5, Nd isotopes with the neutron deficit (A=134-140) and excess (152-158) are 
explored. It is interesting that most of the isotopes are deformed and exhibit prolate shape. 
The cartesian quadrupole moments for the isotopes are given in Fig. 6. This figure also 
shows E1(T=1) centroids and widths. The latter obviously correlate with the magnitude 
of the deformation. The resonance energies generally agree with estimation E = 81A~^^^ 
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FIG. 5: The dipole giant resonance in Nd isotopes with A=134-140 and 152-158, calculated for the 
force SLy6 with the current contribution. The full strength (solid curve) is supplemented by /i = 
and /X = 1 branches (dotted curve) . 

MeV. 

Finally, Fig. 7 demonstrates that the current density contribution to E1(T=1) ends of the 
isotope chain is the same as in the middle of the chain (compare Figs. 4 and 7). So j-impact 
is fully determined by the Skyrme force and does not depend on the shape and neutron 
number of the isotope. Since j-impact is related with the effective masses, this observation 
sends a promising message that results for time-odd densities and effective masses obtained 
for standard nuclei may be relevant for exotic areas as well. 

V. CONCLUSIONS 

The effect of time-odd densities was discussed and examined for the particular case of the 
current density. The contribution of this density to E1(T=1) and E2(T=0) giant resonance 
was explored in the chain of Nd isotopes for a representative set of 8 Skyrme forces. The 
current impact was found to be generally strong and fully determined by the isoscalar and 
isovector parameters Bq and Bi of the Skyrme forces, responsible for the effective masses. 
The results obtained allow to classify the Skyrme forces into 3 groups, depending on the 
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FIG. 6: Energies (upper plot) and widths (middle plot) of the dipole resonance, calculated (SLy6, 
with the current contribution) in Nd isotopes with A=134-158. The bottom plot exhibits the 
calculated quadrupole moments of the isotopes. Dotted line at the upper plot gives the empirical 
estimation E = 81^-^/^ ^eV for E1(T=1) energy. 




FIG. 7: The dipole resonance calculated with the force SLy6 for the cases with (solid curve) and 
without (dotted curve) contribution of the time-odd current in the isotopes ^^^Nd and in ^^^Nd. 



magnitude and sign of Bi. Our study show that time-odd densities represent a principle 
ingredient of Skyrme forces, responsible for the Galilean invariance. They are closely related 
with the effective masses and significantly influence nuclear dynamics. So, time-odd densities 
can be an important factor for further upgrade and unification of Skyrme forces. 

The calculations show that the current impact is fully determined by the Skyrme force and 
does not depend on the shape or neutron number of the isotope. In other words, the impact 
is the same for standard and exotic nuclei. This message hints an interesting possibility to 
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explore the time-odd density effects for standard nuclei (where there are experimental data) 
and then to transfer the results to exotic areas. 
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